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With the new measurements of J/tp and ?/>' -^ yr/c {'yrj'^) from CLEO and BES-III Collaboration, 

we re-investigate the intermediate meson loop (IML) contributions to these radiative decays in 

association with the quark model Ml transitions in an effective Lagrangian approach. It shows that 

^ I' the "unquenched" effects due to the intermediate hadron loops can be better quantified by the new 

D ' data for J/i/> — >■ jrjc. Although the IML contributions are relatively small in J/tp -^ 777c, they play 

^0 , a crucial role in %p' —>■ jrjc (jric)- A prediction for the IML contributions to ?/>(3770) — ^ 777c {"fVc) is 

made. Such "unquenched" effects allow us to reach a coherent description of those three radiative 

transitions, and gain some insights into the underlying dynamics. 
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I. INTRODUCTION 



i-^' There has been a long-standing puzzle on the radiative transition rates of J/ip and tp' — >■ 777c {tt]'^)- 

On the one hand, the nonrelativistic potential model (NR model) including color Coulomb plus linear 

^S| . scalar potential, and spin-spin, spin-orbit interactions, has made great successes in the description of 

^ ' the charmonium spectrum based on the constituent degrees of freedom [15. Nevertheless, a relativised 

C"^ , version by Godfrey and Isgur (GI model) also offers a reasonably good description of the hadron spectra 

■^-Tj ' and transition matrix elements for quarkonia made of either light or heavy 99 [3, |3|. On the other hand, 

both NR and GI model have predicted relatively larger branching ratios for J/ip and ip' — ^ 777c {'^rj'c)- ^^ 

particular, the predicted partial decay width for ip' — >■ 777c was nearly one order of magnitude larger than 

r^ I the experimental data Q . In contrast with the success of the quark model in the description of various 

properties of charmonium spectrum, such discrepancies seem not to be trivial and have initiated a lot of 

theoretical interests in the literature [5l-[l3|. 

Recently the lattice QCD (LQCD) calculations of the charmonium radiative transitions were re- 
ported [IJj ll5|. As shown by Ref. jU], in the "quenched" approximation the magnetic dipole (Ml) 
transition of J/ip -^ 'yrjc was consistent with the new experimental data from CLEO collaboration |16| . 
; ^ ' although one notices that the lattice value does not overlap with the experimental uncertainties. For 

C^ I t/i' — > 777^, the LQCD uncertainties are even larger than the experimental ones, from which one cannot 

conclude that "unquenched" effects would not play a role here. 

In Ref. [12!|, the Ml transition of J/tA — > 7r/c was investigated in the framework of nonrelativistic effec- 
tive field theory. By assuming the ground state charmonium to be a weakly coupled system, the authors 
obtained the radiative decay width T{J/tp — >• 777c) — (1.5 ± 1.0) keV to the correction of 0{v^/ml) with 
a large uncertainty due to high-order corrections. For ip' — >■ 777c(7?7c), the weakly-coupled-cc assumption 
cannot be applied. Nevertheless, the mass of ip' is close to the open DD threshold, which may have 
non-negli gibl e effects on the constituent quark potential. 

In Ref. [l3l, we proposed to consider the intermediate meson loop (IML) corrections as an "unquenched" 
mechanism in the charmonium energy region. Such a mechanism turns out to be important for exclusive 
transitions especially when the mass of the initial state is close to the open channel threshold [1»-^33|. An 
evidential observation of the IML contributions should be the 7/;(3770) non-DD decays. Since 7A(3770) 
is close to the DD threshold, the DD rescatterings into li ght hadrons would be an essential process 
contributing to its non-DD branching ratios [271 . In Ref. [13|, it was shown that the J/tA exclusive 
decays would experience relatively smaller open charm effects than the ip' since the latter is much closer 
to the DD threshold. In another word, the IML would have more important impact on the ip' decays, 
while the J/ip suffers less. 

During the past two years, important progresses have been achieved in the experimental measurements 



of these radiative transitions. The CLEO CoUaboration reported the branching ratios BR{J/ip — > 777c) = 
(1.98 ± 0.09 ± 0.30)% and BR{^' -> 777^) = (4.32 ± 0.16 ± 0.60) x lO'^ ^, while a search for the decay 
of V'' -^ IVc only led to BR{i;' -^ -fr]'^) < 7.6 x 10"^ (H. In fact, the CLEO data for V' ^ IVc have 
been greatly weighted in the PDG2010 averages, i.e. BR{J/ij -^ jr]c) = (1.7 ± 0.4)% and BR{tp' -^ 
IVc) = (3.4 ± 0.5) X 10"'^ [3]. Very recently, BES-III Collaboration also reported the branching ratio 
BR{ip' — >■ 777^) = (4.7 ± 0.9stat ± 3.0sys) x 10^** J4,] as the first measurement of this decay channel. 

The above progresses in both theory and experiment thus prompt us to revisit this problem based on 
the following considerations and improvements of our calculation: i) As shown by the CLEO data [16] and 
LQCD unquenched calculations |14j], there leaves only small corrections from the "unquenched" effects 
in J/tjj — >■ 777c. This will impose strong constraints on our model parameters. As mentioned earlier that 
the IML should have larger impact on the ip' decays, the interest is to investigate how important the 
intermediate D meson loops would be in the ^' channel, ii) In Ref. [l3|, the D*'^ — >■ D^^ coupling was 
fixed by the experimental upper limit for the D*^ total width. This should have overestimated the loop 
contributions involving D*'^D'^j vertex. In this work, we adopt a more realistic coupling value in the 
calculation, iii) It was noticed in Refs. [3^, [Sg that the D*D*{D) loop were rather important. Thus, 
we include all the S'-wave D mesons in the loop calculations to estimate the leading contributions from 
the IML. iv) We do not include the contact terms in this work. Such contributions, though turned out 
to be negligibly small, were considered in Ref. |17| . The contact terms were induced by EM minimal 
substitution at the VVP vertex. We shall discuss later that the contact terms would be eliminated in 
the ELA in order to avoid unphysical contributions, v) The IML transitions also provide a mechanism 
to evade the quark model selection rule at leading order for the Ml transition between a D-wave and 
S'-wave states. Thus, the IML contributions in ■0(3770) should be investigated. 

This paper is organized as below. In Sec. II we present the framework of the effective Lagrangian 
approach (ELA) for the IML. Section III is devoted to numerical results and discussions. A brief summary 
is given in the last section. 

II. EFFECTIVE LAGRANGIAN APPROACH FOR THE IML 

The IML transitions, or known as final state interactions (FSI), have been one of the important non- 
perturbative transition mechanisms in many processes [18|-|32|. In the energy region of charmonium 
masses, with more and more data from Belle, BaBar, CLEO and BES-III, it is widely recognized that in- 
termediate hadron loops may be closely related to a lot of non-perturbative phenomena observed in exper- 
iment, e.g. apparent OZI-rule violations t21j-;32,], sizeable \\o\\-DD decay branching ratios for ■0(3770) [27j , 
and the helicity selection rule violations in charmonium decays J28l - t30| . The IML transitions play a role 
as "unquenching" the simple cc picture in the quark model. It can be easily understood since we know 
that the charm quark is somehow not heavy enough. The failure of the heavy quark approximation will 
then manifest itself in some exclusive transitions, such as the problems investigated in this work. 

Before proceed to the IML formulation, we recall the Ml transition amplitude based on the constituent 
scenario as the following: 

where n{n'), S{S'), L{L'), J{J'), |V'i)(lV'/)) are the initial (final) state main quantum number, spin, 
orbital angular momentum, total angular momentum and spatial wavefunctions, respectively. E^ and 
Ef denote the final state photon and meson energy, respectively, while Mi is the initial cc meson mass. 
The above equation can be regarded as the "quenched" contributions. The GI model Ml radiative rates 
do not incorporate the phase factor Ef/Mi, while include a recoil factor ja{kr/2). In this work, we will 
quote the results of Ref. Q . 

In the VVP transition, all mechanisms that contribute to this transition will appear as corrections to 
a single anti-symmetric tensor coupling. We derive the effective VjP couplings via the Ml transition of 
Eq. ([1]) as follows, 

Mf^{Ml) = gvjpeo.p^,pfe^p';e; , (2) 



where pi and p^ are four- vector momentum of the initial meson and final state photon, respectively, and 
£i and £-y are the corresponding polarization vectors. 

The IML transitions can be schematically illustrated by the triangle (Fig. [1]) and contact processes 
(Fig. [5]). The coupling between an 5*- wave charmonium and charmed mesons is given by the effective 
Lagrangian based on heavy quark symmetry |37l . |38| , 

C2 = ^g2Tr[R,sH2^l''^^.Hl^] + H.c, (3) 

where the iS-wave vector and pseudoscalar charmonium states are expressed as 

Rc-c = (^ii^) in, ^c75) (^) . (4) 

The charmed and anti-charmed meson triplet are 

Hu = (^)[2?*%-A75], (5) 

H2. = [A*^p-^.75](^), (6) 

where V and V* are pseudoscalar ((_D*',_D+, !?+)) and vector charmed mesons ((I?*°, I?*+,-D*+)), re- 
spectively. 

Explicitly, the Lagrangian for the S-wave charmonium (J/ip, ip' and rjc) couplings to D and D* becomes 

+g^^D*DD*^^{d,,r^,D - 7^,d,,D) + igr,^D'D'e^''"^^d^DlD*^dpT^, . (7) 

The effective Lagrangians for the electromagnetic (EM) interaction vertices oi^DD, ^D*D* and 'yDD* 
are [39| 

CoD-yix) = ieA^{x)D-{x)^''D+{x)+ieAf,D-{x)^''Dt{x) (8) 

Co^n-^ix) - teA^ix) L'^^ Dl-^f^D^+ix) + g^^ D:-{x)d'^D;+{x) ^ g'^'^ df'D:-{x)D;+{x) 

+ 7eA,{x) |g"^ Dl^'^^Dl+ix) + g^^ D*,- {x)d" 0*+ {x) - g'^" d^ DI-{x)D:+{x) \ , (9) 



Cd'D^{x) = -e^''"'PF^^{x)igD^-D+^Dl;-pix)D+ix)+gD'Ooo^D*„''p{x)D^{x) 
+ 9DrDt->D:-0ix)Dtix)^+ H.c, 



(10) 



where A d f,B = Adf,B- d^AB, F^,y = d^Ay - d^A^, and M^^, = d^M^, - dyM^,. The EM interaction for 
neutral D meson [D^D^^ and I?*°Z?*°7) do not exist. The coupling constants appearing in the effective 
Lagrangians will be determined later. 

The loop transition amplitudes in Fig. [T] can be expressed in a general form in the ELA as follows: 

Af MV2 V- T1T2T3 ^ 

J (27r)^ ^^^ 010203 

D' pol. 

where Ti {i — 1,2,3) are the vertex functions; ai = pj — mf (i — 1,2,3) are the denominators of the 
intermediate meson propagators. We adopt the typical dipole form factor in the calculation, i.e. 

■^(-2,^^)^(^^)', (12) 
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FIG. 1: Schematic picture for the decay of J/ip -^ -yrjc via (a) D+D~{D*+), (b) D+D*'{D+), (c) D+D*-{D*+), 
(d) D*+D~{D*+), (e) D*+D*~{D+), (f) D*+D*"(D*+) intermediate charmed meson loops. Similar diagrams 
are for strange charmed mesons. For neutral charmed mesons, only (a), (c) and (e) can contribute. Similar 
pictures occur in tp' — > 77^0 and 777^. 



where A = 1112 + aAgcD and the QCD energy scale, Aqcd = 220 MeV. This form factor will take care 
of the non-local effects of the vertex functions and kill the loop divergence in the integrals. The value of 
parameter a is commonly at the order of unity. 

In Fig. [TJ the triangle diagrams for charged intermediate meson loops are illustrated. For the neutral 
ones, only those corresponding to diagrams (a), (c) and (e) can contribute. Diagrams (b), (d) and (f) 
have no contributions due to the vanishing D^D^^ and D^*D*^^ couplings. Note that it is easy to check 
that all the diagrams in Fig. [T] satisfy gauge invariance individually. The explicit transition amplitudes 



for those triangle loops are given as follows 

(2tt 



■^T = (*^) / 77;::^[9^DDe''^{pip- PZp)\[-gy^^D'D{Pfa +P3.a)][-egD*D'iealif,„p'^eZJ^p':^ 



,^^^ 



- -2 2 2 2 2 2^i'^2,P2) , (13) 

pj -mf p^- m^ p^ - ml 

7^-^[5»A-D-I)ea/3M^P?^^b3 -Pl)]h5»fc-D--D(P/p +P2p)][ee*''(pia +P2^)] 

^ p\-m\pl-ml pl-ml -^("^2,^2), (i4j 

/74 

.. ^ »(-g"^ +P2P2/^i) »(-g^^ +P3P3/"^i) ^. 2^ ., .^ 

2 13 l2 112 l2 112 -^l"^2,P27 J U^j 



■^T^ = i^^) 77rT^h94'D*DSaPp^P^e'^^{P3-Pi)][-gvaD'-D{pfe+P3e)] 



Pi - ™l Pi - ^i P5 

(27r) 

[ee7(~5ffT(Plp + P2p) + gprPla - gpaP2T)] 

ii-gf'^+p'^pl/ml) i{-g''^ ^p%pllml) 



- 2 2 2 2 2 2-^("^2,P2) > (16) 

p^ — m^ P2 ~ "^2 P3 ~ "^^3 

/74 
7^-^£j^[ff0£)*D-((Plp -P3p)gali + .g/3MP3a " .ga^Pl/s)] [-g,fc_D-D(P/,. + P3i.)] 

X f-ean.n £ . .p.^^.Ci ^('g"^ + P?K/m^) » »(-/-+ pgpgMi) -..„^^ 2. .^7. 

X [ e-gu'D-iS-pai-rP^e^ Px\ „2 „ ™2 „2 „ ™2 „2_™2 -^V™2,P2j; U'j 



/74 
7^:^evl[5V£'*£"((PlM --P3M)5a/3 + gUpPZa ~ gapPip)][gr,,D*D'epcrirP2P}] 

r *(9/ / N ^Ji-g°"^ +PlPl/mf) 

X [-ee^ (.90k(P20 +Pie) + 5fc0Pi0 - 500P2k)J 2 2 

Pi "^1 

"^ ^f^^;^! ^F^^^l -^("^2,^2) • (18) 

The contact diagrams of Figs.[2ja), (b), (c) and (d) arise from gauging the strong J/'ipD*D and t]cD*D 
J/ipD*D* and r]cD*D* interaction Lagrangians by the minimal substitution d^ ^ d^ + ieA'^. One can 
also easily check that gauge invariance is guaranteed for the contact diagrams in Fig. [2j However, these 
contact diagrams can be neglected in the calculation based on the following detailed examinations. The 
processes in Figs. [2] can be classified into two categories. The first one includes diagrams (a) and (d), 
where the contact vertices induced by the EM minimal substitution violate gauge invariance. We discard 
these contributions based on the empirical argument that an initial massive vector meson decaying into 
jVP via the contact interaction actually violates gauge invariance, thus is forbidden. The second category 
includes diagrams (b) and (c) in Figs. [2] Although these contact vertices keep gauge invariant, they do 
not have contributions due to the vanishing loop integrals as shown in Ref. jlTJ ]. In brief, we argue that 
the contact diagrams would not contribute to the transition matrix elements in the ELA. Therefore, we 
can concentrate on the triangle diagrams in the following calculations. 

In principle, we should include all the possible triangle meson loops in the calculation. In reality, 
the breakdown of the local quark-hadron duality allows us to pick up the leading contributions as a 
reasonable approximation [40, |4l|. Also intermediate states involving flavor changes turn out to be 
strongly suppressed. One reason is because of the large virtualities involved in the light meson loops. 
The other is because of the OZI rule suppressions. So we will only consider the charmed meson loops as 
the leading contributions in this work. It should be noted that the IML transitions can naturally evade 
the quark model selection rule as a dynamical mechanism. 
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FIG. 2: The contact diagrams considered in J/ip — >■ 777c. Similar diagrams are also considered in 1/)' -> ^ric{ri'i.). 

III. NUMERICAL RESULTS 

To proceed the numerical results, we first clarify the following points: 

(i) In the heavy quark limit, the coup lings for charmonium-charmed mesons in Eq. (I7|) can be related 
to the parameter 172 defined in Ref. [37|, 



g-ipvv = 2g2y/m^mv, g-iiivv = -~2g2^/m^rnv' , g^jivv ='2g2\ , (19) 

where 52 = y/fntii / {'^rnxi fip) and the /^ — 405 MeV is the J/ip decay constant. 

(ii) The ratios of the couplings constants g^'OD and g-ip'O'D* to gj/tpoD and gj/tpo'D' are fixed as 

g^^D^ = 92^21^^ = 0.9 . (20) 

gj/i)DD gj/tpD'-D" 

This ratio has uncertainties as adopted in the literature [30|. A recent study of the e+e" — >■ DD cross 
section lineshape [35| suggests that g^p'OD — gj/tpDD, and we will discuss later that to coherently account 
for the partial widths of J/^ and -ip' — )■ 777c (jr]'^) would impose a strong constraint on the ratio of 

gtp'Do/gj/ipDD- 

(iii) The radiative couplings for D* — > D^ can be determined by the partial widths from experimental 
measurement, i.e. 

r{D*^Dj)^^gl,^^\P,\^ . (21) 

The partial width r(D*+ -^ D+^) has been precisely measured, i.e. T{D*+ -^ D+^) = 1.54 keV [3]. 
This allows us to extract gD*+D+-y = —0.5 GeV~^. For £)*° — >• D°7, the branching ratio is measured 
by experiment [J|. However, the total width of D*° has not been well determined. This will bring some 
uncertainties to the estimate of the radiative coupling. However, this quantity can be related to the 
g£)*+D+^ in the constituent quark model, which gives gD*'>DOj = +2.0 GeV~ . This is also a value 
obtained in different approaches. For the coupling go^Ds-y^ the value (—0.3 ± 0.1) GeV~^ from QCD 
sum rules (QSR) [42] is adopted. We note that their relative signs are consistent with each other in the 
framework of LQCD [43], QSR [43], and constituent quark model. 

(iv) Another two undetermined parameters in our model are the cut-off parameter a in the form 
factor and the relative phase S between the "quenched" (i.e. quark model Ml transition amplitude) 
and "unquenched" (i.e. IML transition) amplitude. Taking the advantage of the anti-symmetric tensor 
coupling for VVP, we can always parametrize the total amplitude as 

Mf^ = [gvjp + ~gtr^e'^]eo,p,,,p'^e^^p'^e\ , (22) 
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TABLE I: Radiative partial decay widths given by different processes are listed: rj^f and Vfj^ are the Ml 
transitions in the NR and GI model, respectively ,3]. Viml are contributions from the IML transitions (Fig. [T]), 
and Fan are coherent results including the Ml transition amplitude of the GI model and IML transitions. The 
experimental data are from PDG j^l, BES-III [3^], and CLEO [IGj. The LQCD results are also included as a 
reference llal. 
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FIG. 3: Partial widths of i/)' — >■ 777^ in terms of the mass of 77^ with q = 0.98. The dashed line is the exclusive 
IML contributions, and the dotted line is the exclusive GI model result. The solid line is the combined result of 
both. The datum point is from the BES-III Collaboration [34 ]. 



where gv^p is fixed to be a positive real number by the "quenched" quark model Ml transition amplitude, 
and gtri is extracted from the triangle loop diagrams. We simply take (5 = or tt in the calculation, since 
the decay threshold of the intermediate mesons are above the initial meson ( J/t/) and '0') masses, the 
absorptive part of the loop integrals does not contribute as a consequence. In fact, because the "quenched" 
quark model Ml transition amplitude has overestimated the experimental data, it determines (5 = tt in 
the calculation. 

Note that in the "quenched" quark model scenario the spin-flipping Ml transition amplitude for J/f/' -^ 
^r\c shares the same sign as that for ?/;' — >■ 777^ since the latter pair states are just the corresponding radial 
excitations of the former ones. In this sense, the same phase angle (5 = vr for -0' -^ JVc seems to be 
reasonable. It is interesting to observe that this phase relation is also respected in ip' — >■ 777c- 

With the destructive phase S = tt, the parameter a = 0.98 ± 0.27 is determined by combining the GI 
model result with the IML to reproduce the experimental partial width TexpiJ/tp ~^ IVc) = (1.58 ±0.37) 
keV Q. The range of a is given by the experimental error bars. We then apply the same set of a and 6 
to predict the partial width of ^p' — ?> 777c(7??c)> ^^'^ "unquenched" effects in ^(3770) — > 777c and 777^. 



In Tab. U the calculated Ml transition branching ratios are listed to compare with the GI and NR 
quark model results. The partial width of J/"!/; — > 777c is an input to fix a. The LQCD calculations are 
also included as a reference. In comparison with the previous estimate in Ref. |17| . the requirement of 
smaller "unquenched" effects leads to smaller IML contributions in all these three decay channels. 

For -0' — > 7?7c {jT]'^)^ we learn the following points: (i) The "unquenched" contributions from the 
IML still play a role of destructive interferences with the GI amplitude to bring down the Ml transition 
amplitudes. The calculated partial decay widths are consistent with the BES-III preliminary result of 
T{tp' — >■ 777^) [3JI, and CLEO measurement of T{ip' -^ ^rjc) [i^|. (ii) One notices that the uncertainty 
with the form factor parameter a extracted in J/ip — ^ 777c still causes large uncertainties in the estimate 
of the IML contributions in the ip' decays. This shows a sensitivity correlation of the IML contributions 
for the n^Si — )■ n'^So Ml transitions. One also notices that the partial width r(^/'' — >■ 777^) has relatively 
smaller uncertainties than T(ip' — > 7?7c). This is self-consistent since the former transition does not violate 
the selection rule of Eq. ([T]) and the "quenched" quark model leading order transition is still dominant. 

For 'ijj{3770) — >■ -f7]c (jVc)^ if V''(3770) is a pure D-wave state, the Ml transition will be forbidden by the 
selection rule of Eq. ([T|). However, due to the nonvanishing photon energy in the decay, higher multipoles 
beyond the leading one would contribute. In a harmonic oscillator basis, the nonvanishing transition 
amplitude of ^^(3770) — ^ IVci'n'c) i^ the same order as that of V'' ^ IVc (Vc)- Since a quantitative estimate 
of the quark model amplitude will depend on the details of model constructions, we only concentrate 
here the IML mechanism that present the "unquenched" contributions. As listed in Tab. HI the IML 
transitions predict r(V'(3770) -^ 777^) = (17.14J:?2:o3) keV and r(V'(3770) -^ jr]'J = {1.82tlil) keV, 
which are in a reasonable order of magnitude, although uncertainties appear to be significant. This 
is similar to ip' — >■ 777c, where sensitivity of the partial widths to the range of a values is obvious. 
Interestingly, it shows that the IML contributions are the same order as the LQCD results [1^1 • This 
implies that interferences between the "quenched" and "unquenched" amplitudes should be important 
for the ip{3770) radiative decays. As a consequence, the radiative transition of V'(3770) could become 
either abnormally strong or significantly small in comparison with potential quark model expectations. 
Experimental measurement of these radiative transitions would be helpful for providing further constraint 
on the IML contributions. 

It is interesting to note that the present experimental data for J/ijj and -0' -^ ^rjc (7?7c) would tightly 
stretch the parameter space for the form factor parameter a and coupling g^'OD- Since these processes 
should share the same form factor parameter a, the main parameter difference is the coupling between 
J/tp{ip') to £)(*)!)(*) for which the analysis of Ref. [S^ suggests that a smaller value for gt/j'DD should be 
applied. 

For ip' — > 77/^, the mass of 77^ may also cause uncertainties to the extracted IML contribution. With 
the fixed a = 0.98, we investigate the sensitivities of the IML contributions to the 77^ mass. In Fig. [3J 
we plot the exclusive partial width T^ip' — >■ 777^) from the IML in terms of the mass of 77^ within a range 
of S.GSS'-^ 3.641 GeV |4]. It shows that within the PDG mass range, the IML contributions are rather 
stable. The partial width decreases in term of the increasing r?z,,' due to the decreasing phase space in 
the decay transition. 

IV. SUMMARY 

We revisited the hadronic meson loop contributions to the J/ijj and ip' radiative decays into 777c or 
777^ in the ELA. In the framework of an improved effective Lagrangian approach, the IML transitions 
provide "unquenched" corrections to the leading couplings extracted from potential quark models due 
to the unique Lorentz structure of VVP interactions. In comparison with the NR and GI model, the 
IML contributions intend to cancel the quark model "quenched" amplitudes. Apart from those more 
elaborate treatments for the meson loop calculations, we have applied an experimental constraint on the 
strong coupling of g^'DD based on the analysis of the cross section lineshape of e+e~ — )• DD [S^, where 
a relatively smaller value of g^p'OD was favored. It is interesting to see that the IML effects in J /ij} — >■ 777c 
are much smaller than that in ■0' — > 77/c iiv'c)- Note that the pure Ml contribution in ip' — > 777c is about 
one order of magnitude larger than the experimental data, such a significant discrepancy implies the 
necessity of "unquenching" the quark model scenario especially when the transitions are close to open 



thresholds. 

For -0(3770) -> jrjc {jri'^), we predict quite significant_corrections from the IML, which are the same 
order of magnitude as the LQCD "quenched" resuh pjj. This is an interesting issue related to the 
■0(3770) non-D5 decays. The BES-III Collaboration recently scanned over the ?/'(3770) mass region. It 
will be possible to measure the radiative decays of "0(3770) — > ^rjc i^i]'^), and further clarify the role 
played by the IML effects as an "unquenched" mechanism for the cc quark model scenario. 
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